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1.  INTRODUCTION 


The  develq)ment  of  residual  stress  In  composite  laminates  is  strongly  influenced  by  processing  history. 
Residual  stress  can  have  a  signiflcant  effect  on  the  mechanical  performance  of  composite  structures  by 
inducing  warpage  or  initiating  matrix  cracks  and  delaminations  (Kau  and  Petrusha  198S;  Kays  1985). 
Processing  concerns  associated  with  thermosetting  composites  become  increasingly  important  for 
components  of  appreciable  thickness.  The  most  familiar  problem  is  an  increase  in  temperature  resulting 
from  the  resin  exothennic  chemical  reaction  that  may  raise  internal  temperatures  to  levels  inducing 
material  degradation.  A  secottd  concern  is  the  complex  gradients  in  temperature  and  degree  of  cure 
accentuated  by  increased  thickness  (Bogetti  and  Gillespie  1991). 

Traditional  residual  stress  analysis  in  thermosetting  composite  laminates  is  based  on  thermal  expansion 
mismatch  between  adjacent  plies,  a  uniform  temperature  drop  from  the  cl'?  temperature  to  ambient 
cortditions.  and  no  stress  development  prior  to  comi^etion  of  the  curing  process  (Griffln  1983;  Hahn  and 
Pagano  1975,  1976;  Stango  and  Wang  1984).  This  approach  may  not  be  adequate  for  predicting  the 
process'induced  stresses  in  thick-section  thermosetting  composite  laminates.  The  authors  have  shown  that 
spatial  solidification  (tesulting  from  complex  temperature  and  degroe-of-cure  gradients)  introduce  several 
important  mechanisms  not  accounted  for  in  the  traditional  stress  analysis  methodology  (Bogetti  and 
GiUespic  1989,  1990a.  l99Qb). 

The  mechanisms  governing  stress  development  in  thick-section  composites  ire  quite  similar  to  those 
encountered  in  the  manufacture  of  tempeo’d  glass  (Lee.  Rogers,  and  Woo  196.';  Tackels  and  Crochet 
1973;  Crochet  et  al.  1974;  Crochet  and  Denayer  19<j0).  iuch  stresses  result  from  tlie  interactions  between 
spatially  varying  thermal  contractions,  viscoelastic  material  response,  and  severe  temperature  gradients 
which  develop  during  the  quenching  process.  Early  treatment  of  stress  development  due  to  viscoelastic 
material  response  in  the  presence  of  thermal  gradients  was  based  on  the  time-temperature  superposition 
principle  of  theologically  simple  material  behavior  (Muki  and  Sternberg  1961).  This  approach  has  been 
applied  to  study  the  macroscopic  in-plane  residual  stress  development  during  the  quenching  of  a  folly 
cured  epoxy  plate  (Maneschy  et  al.  1986)  and,  more  recently,  to  the  quenching  of  thermoplastic  matrix 
composites  (CJhapman  et  al.  1988). 

As  a  thermosetting  resin  cures,  the  effective  mechanical  properties  of  the  composite  vary  from  a 
viscous  liquid  (negligible  stiffness),  in  its  uncured  state,  to  a  viscoelastic  or  clastic  solid  (high  stiffness). 
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in  its  fiilly  cured  state.  During  cure  the  resin  undergoes  significant  increases  in  stiffness  and  volumetric 
shrinkage  associated  with  the  cross-linked  polymerization  reaction.  The  mechanical  properties  of  the  resin 
phase  during  cure  are  governed  by  competing  mechanisms  between  chemical  kinetic  hardening  and 
viscoelastic  relaxation  phenomena.  Once  fully  cured,  tire  resin  exhibits  traditional  viscoelastic  behavior 
at  elevated  temperatures  and  appioaches  elastic  behavior  at  lower  temperatures.  Additionally,  thermal 
expansion  of  the  fully  cured  resin  is  the  only  mechanism  contributing  to  changes  in  specific  volume. 

Complex  temperature  and  degree-of-cure  gradients  which  develop  in  thick-section  thennosetting 
laminates  induce  spatially  varying  material  response  associated  with  resin  stiffening  and  shrinkage  as  the 
chemical  reaction  progresses.  Birefringence  patterns  have  shown  experimentally  that  temperature  and 
degree-of-cure  gradients  have  a  significant  effea  on  the  development  of  residual  stresses  in  thick 
thermosetting  castings  (Pusatcioglu  et  al.  1980).  In  a  series  of  papers,  Levitsky  and  Shaffer  (1974, 1975) 
and  Shaffer  and  Levitsky  (1974)  have  studied  the  developmem  of  residual  stresses,  including  temperature 
gradients  and  spatially  varying  "chemical  hardening”  effecu  on  mechanical  properties  in  isotropic 
thermosetting  materials  exhibiting  no  chemical  shrinkage.  Chemical  s)?tinkage,  however,  represents  a 
potentially  significant  change  in  specific  volume  (superimposed  on  changes  in  specific  volume  due  to 
temperanue)  that  exists  in  theimoset  composites  during  the  curing  process. 

Recently,  Bogetti  and  Gillespie  (1989,  1990a,  1990b)  developed  a  model  to  predict  process-induced 
stress  in  thick  thennoset  laminates  employing  a  cure  simulatipn  analysis  based  on  an  incremental  transient 
finite  difference  formulation  that  accounts  for  thermal  and  chemical  interactions  during  processing. 
Material  models  were  proposed  to  describe  the  resin  modulus  and  volumetric  (chemical)  shrinkage  as  a 
function  of  degree  of  cure.  A  micromechanics  model  for  continuous  fiber  composites  was  used  to  evaluate 
the  instantaneous  spatially  varying  mechanical  properties,  thermal  expansion,  and  chemical  shrinkage 
strains  through  the  laminate  thickness  as  a  function  of  temperature  and  degree  of  cure.  Process-induced 
stress  and  deformation  predictions  are  based  on  an  incremental  laminate  theory  that  includes  temperanire 
gradients,  spatially  varying  cure-dependent  mechanical  properties,  thermal  expansion,  and  chemical 
shrinkage  strains.  Details  of  the  cure  simulation  and  process-induced  stress  models  are  presented 
elsewhere  (Bogetti  and  Gillespie  1989,  1990a,  1990b,  1991). 

The  models  developed  demonstrate  the  complex  interactions  that  exist  between  processing  conditions 
(thermal  history),  cure  kinetics,  resin  shrinkage,  laminate  thickness,  and  stress  development  in  thick 
thermosetting  composite  laniiruktes.  Previous  investigations  have  focused  on  the  influence  of  thermal 
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history  and  laminate  diickness  on  stress  development  in  thick  glass/jpolyester  laminates.  In  this  work,  the 
ptocess'induced  stress  model  is  onpioyed  in  parametric  studies  aimed  at  assessing  the  influence  of  resin 
chemical  shrinkage  and  thermal  expansion  on  stress  development  in  thick-section  glass^lyester  and 
gra[Mte/epoxy  laminates.  Model  predictions  are  also  compared  with  previously  published  experimental 
data  for  cure-dependent  modulus  and  dimensionless  laminate  curvature  of  a  graphite/epoxy  (AS4/3S01-6) 
cross-ply  laminate  (Hahn  and  Kim  1990). 

2.  ANALYSIS 

2.1  Cure  Simulation  and  Incremental  Stress  Analysis.  The  process-induced  stress  model  employed 
in  this  investigation  consists  of  three  key  submodels:  (1)  a  cure  simulation  analysis,  (2)  material  models 
to  describe  the  resin  behavior  during  cure,  and  (3)  an  incremental  laminate  plate  theory  stress  analysis. 
The  theoretical  development  and  integration  of  these  submodels  is  not  presented  here,  as  detaiis  are 
documented  elsewhere  (Bogetti  and  Gillespie  1989. 1990a,  1990b,  1991).  Since  this  investigation  focuses 
on  the  effects  of  resin  behavior  during  cure  (chemical  hardening  and  shrinkage)  on  process-induced  stress 
development,  only  these  material  models  are  reviewed  in  this  report 

2.2  Material  Models.  The  incremental  stress  analysis  requires  as  input  the  mechanical  properties  and 
macroscopic  strains  of  the  composite,  undergoing  spatial  chemical  hardening  and  shrinkage  in  a  transient 
thermal  environment.  The  effective  composite  mechanical  properties  and  strains  are  based  on  a 
unidirectional  composite  micromechanics  model  The  constituent  fiber  properties  (constant),  the  resin 
properties  and  chemical  shrinkage  (cure  dependent),  and  the  fiber  volume  fraction  are  input.  Two  material 
models  are  proposed  to  describe  the  mechanical  properties  and  volumetric  shrinkage  of  the  resin  during 
cure.  Changes  in  the  resin  properties  directly  influence  the  mechanical  properties  in  the  composite,  and 
chemical  shrinkage  represents  a  significant  source  of  internal  loading  during  cure. 

The  first  material  model  defines  the  instantaneous  resin  modulus,  E„.  explicitly  in  terms  of  the  degree 
of  cure: 


-  (1  -  a)  V  +  a£,oo  ♦  yad  -  a)(E,oo  -  e^°)  . 


(1) 
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where  and  are  die  assumed  liiUy  uncured  and  fully  cured  resin  moduli,  respectively.  The  tenn 
y,  which  must  satisfy  (- 1  <  y  <  1),  is  introduced  to  quantify  the  competing  mechanisms  of  stress 
relaxation  and  chemical  hardening  (Dillman  and  Seferis  1987).  Increasing  y  physically  corresponds  to  a 
more  ra^d  increase  in  modulus  at  lower  degree  of  cure  before  asymptotically  approaching  the  fully  cured 
modulus.  A  value  of  y  s  0  was  used  to  generate  all  the  results  presented  in  this  woik.  The  instantaneous 
resin  shear  modulus  (also  required  input)  is  based  on  the  isotropic  material  relation  with  v.  assumed 
constant: 


■ 


(2) 


The  inherent  limitation  of  this  material  model  is  diat  stress  relaxation  cannot  be  modeled  explicitly 
(i.e.,  modulus  is  monotonically  increasing  with  cure).  Stress  predictions  resulting  from  this  model, 
therefore,  provide  an  upper  bound  estimate  on  process-induced  stresses  to  be  expected  during  cure.  It  is 
anticipated  that  this  model  would  be  most  relevant  to  material  systems  undergoing  rapid  cure. 

The  second  material  model  describes  the  volumetric  chemical  shrinkage  of  the  resin  during  cure.  No 
chemical  shrinkage  occurs  after  the  cure  is  complete.  The  resin  chemical  shrinkage  induces  significant 
macroscopic  strains  in  the  composite,  representing  an  import^  source  of  imemal  loading  in  thick-section 
laminates  in  addition  to  the  traditionally  recognized  thermal  expansion  strains. 

Neglecting  higher  order  terms,  the  isotropic  resin  chemical  shrinkage  strain  of  a  unit  volume  elemern 
of  resin,  Ae,,  resulting  from  an  incremental  volume  resin  shrinkage,  Av^,  is  given  by 


-1  . 


(3) 


The  incremental  volume  resin  shrinkage  is  based  on  an  incremental  change  in  degree  of  cure.  Aa,  and  the 
total  volumetric  shrinkage  of  the  completely  cured  resin,  V  from  the  expression 

Av,  ■  AaV*^*  .  (4) 
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Incfcmental  effective  longitudinal  and  transverse  chemical  shrinkage  strains  in  the  composite, 
and  respectively,  are  computed  over  each  time  increment  in  the  cure-simulation.  These  strains 

are  based  on  a  micromechanics  model  (Hill  1965;  Whitney  and  McCullough  1980)  utilizing  constant  fiber 
properties,  instantaneous  cure-dependent  properties  of  the  resin,  the  resin  chemical  shrinkage  strain 
increment,  Ae,,  a  zero  chemical  shrinkage  strain  in  the  fiber,  arkl  the  fiber  volume  fraction. 

Incremental  effective  thermal  expansion  strains  are  also  computed  over  each  time  increment  during 
the  simulation.  These  strains  are  based  on  temperature  increnrents  between  time  steps  in  the  simulation, 
AT.  and  the  instantaneous  effective  longituditud  and  transverse  thermal  expansion  coefficients  of  the 
composite,  artd  07-,  respectively.  The  incremental  longitudinal  and  transverse  strain  increments  in  each 
ply  are  defined  as: 


Ae^t  -  C4.Ar 

AeS  •  ct^r  .  (5) 


The  expansion  coefficients,  and  are  based  on  the  micromechanics  model  utilizing  constant  fiber 
properties,  cure  dependent  resin  properties,  constant  thermal  exparuion  coefficients  of  both  the  fiber  and 
resin,  as  well  as  fiber  volume  firactiotL 

The  total  stress-free  macroscopic  processing-induced  ply  strains  are  computed  from  the  superposition 
of  the  thennal  and  chemical  strain  increments  during  the  curing  process.  Gradients  in  the  material 
response  to  chemical  shrinkage  and  thermal  expansion  strains  during  cure  hardening  represem  significant 
sources  of  internal  loading  accounted  for  in  the  stress  model  employed  in  this  investigation.  Details  of 
the  material  models  presented  here  and  the  process-induced  stress  analysis  methodology  are  presented 
elsewhere  (Bogetti  and  Gillespie  1989,  1990a,  1990b.  1991). 

3.  RESULTS 

Required  input  into  the  stress  model  is  rather  extensive;  consequently,  an  input  parameter  summary 
is  presented  first.  Model  predictions  for  resin  modulus  development  and  volumetric  shrinkage  of  the 
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composite,  associated  with  resin  ctossUnUng  during  isothermal  cure,  are  presented  for  glass/);x)lyester  and 
graphite/epoxy.  Cure  simulations  are  tun  at  several  temperatures  for  each  system  to  examine  the  effects 
of  crue  kinedcs  with  regard  to  the  effects  of  signiiicam  thiough-the-thickness  temperature  gradient^  The 
iiUluence  of  chemical  volumetric  shrinkage  on  the  residual  stress  pioflks  in  both  [^ass^wlyester  and 
graphite/epoxy  laminates  is  examined.  Finally,  correlation  of  model  predictions  of  transverse  modulus 
development  in  unldlrecdonal  laminates  and  of  curvature  development  in  an  unsymmetrical  [0/90] 
graphite/epoxy  laminate  is  made  with  previously  pubtished  experimental  results. 

3.1  Input  Parameter  Summaiv.  Irqmt  parameters  used  in  this  investigation  are  summarized  in  this 
section.  The  thermal  properties  of  the  composite  systems  are  presented  in  Table  1.  The  complete 
descripdon  of  the  cure  kinetics  for  the  composite  iiKludes  the  total  heat  of  reaction  and  a  description  of 
the  rate  of  reaction  as  a  function  of  temperature  and  degree  of  cure.  Reaction  rate  expressions  for  the 
glass/polyester  ard  graphite/epoxy  material  systems  are  different  in  form  due  to  the  inherem  differences 
in  the  overall  order  of  the  reaction  kinetics. 


Table  1.  Thermal  Properties  for  Olasa/Polyester  and  Oraphite/Epoxy 


Material 

p[kg/ta’] 

CpIkJ/w'C) 

k,(kW/(m*0) 

Glass/Polyester 

Oraphite/Epoxy 

1.89  X  10* 

1,52  X  10* 

1.26 

9.42  X  la* 

2.16  X  10^ 

4.46  X  10^ 

The  glass^lyester  material  system  consists  of  CYCOM  41Q2  polyester  resin,  manufactured  by  the 
American  Cyanamid  Corporation,  and  is  reinforced  with  E'glass  hbers  (54%  by  volume).  The  reaction 
rate  expression  for  the  glass^Iyester  system  is  second-order  overall  in  the  sense  that  (m^  ♦  n,  •  2)  (Adams 
1988); 


-  A.  exp  (-AE,//r/)cr'(  1  -  a)"*)  .  (6) 

The  constant  is  the  universal  gas  constant  and  7  is  absolute  temperature.  The  exponents  and  n,. 
the  preexponential  coefficient,  A^,  the  activation  energy,  and  the  total  heat  of  reaaion  are  listed  in 
Table  2. 
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Table  2.  Cure  Kinetic  Parameters  for  Glass/Polyester  and  Gr^)hite/Epoxy 


Glass/Polyester 

1 

0.524 

1 

1.476 

1  AJmin,-!) 

3.7  X  10“ 

1  AEJJAnol] 

1.67  X  10* 

I  H,[kJ/kgl 

77.5 

Graphiu 

iiEpoxy 

A,[mia-1] 

2.102  X  l(f 

A2[mitL-l] 

-2.014  X  lO* 

Ajimia-l] 

1.960  X  10* 

AE,[J/mol] 

8.07  X  10* 

AEjtJAnoI] 

7.78  X  10^ 

AEjtJ/mol] 

5.66  X  10^ 

H,(kJ/kg] 

198.9 

The  graphite/epoxy  material  system  contains  Hercules  Corporation’s  3501-6  resin,  reinforced  with 
unidirectional  AS4  graphite  fibers  (67%  by  volume).  The  reaction  rate  expression  for  the  giaphite/epoxy 
system  follows  a  markedly  different  form  (Bogetti  and  Gillespie  1991): 


_  -  ik,  *k,a){l  -0(0.47  -O  M  (aiO.3) 
dt 

and 


^•*,(1-0  /or  (a  >0.3)  (7) 


where  kj,  arxl  ikj  are  defined  by  the  Arrhenius  rate  expressions: 


ik,  ■  Afixp{-6E^IRT) 

ikj  ■  A,exp(-A£^^ 

k^  -  A,exp(-A£^/?r)  .  (8) 
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The  preexponendal  coefficients  >4 p  A,,  and  A,,  the  activation  energies,  AC,.  ACj,  and  ACj,  and  the  total 
heat  of  reaction  for  the  graphite/epoxy  composite  are  also  summarized  in  Table  2. 

The  fiber  and  matrix  constituent  medianical  properties  employed  in  all  the  parametric  studies  are 
summarized  in  Tables  3  and  4. 

The  results  generated  in  this  investigation  arc  sensitive  to  both  the  time  step  size  and  the  grid  spacing 
used  in  the  simulation.  It  is  noted  here  that,  for  the  range  of  processing  variables  and  materials  specific 
to  this  investigation,  time  step  size  and  grid  spacing  incremetus  were  selected  sufficiently  small  to  yield 
converge  solutions  for  all  the  example  problems  presented.  Specifically,  a  time  step  size  of  2.0  s  and  a 
grid  spacing  increment  of  0.00127  m  were  used  to  generate  all  the  results  presented  in  this  investigation. 

3.2  Chemical  Kinetic  Effects  on  Materia*  Models.  Figure  1  shows  the  development  of  dimensionless 
resin  modulus  vs.  time  during  isothermal  cure  for  both  the  glass^lyester  and  graphite/epoxy  systems. 
Isothermal  cure  isolates  the  behavior  of  the  system  induced  by  chemical  kiitetics  from  that  induced  by 
transient  thermal  effects.  Temperatures  at  which  the  simulations  were  run  correspond  to  a  range  of  the 
respective  manufacturer’s  suggested  core  temperatures.  Upon  comparison  of  the  curves,  it  is  apparent  that 
the  cure  of  the  polyester  resin  is  much  more  violent  than  that  of  the  epoxy.  More  importantly,  we  see  that 
the  polyester  system,  due  to  the  extremely  lapid  nature  of  its  cure  behavior,  is  inherently  much  more 
sensitive  to  temperature  variations  than  the  epoxy  system*.  Figure  2  illustrates  the  development  of 
chemical  volumetric  shrinkage  resulting  from  the  cross-link  polymerization  reaction  for  both  systems. 
Again,  the  rapid  nature  of  the  polyester  cure  with  respect  to  the  epoxy  resin  is  noted  (shrinkages  have  been 
nondimensionalized).  We  see  that  small  temperature  variations  in  the  glasi^lyester  system  may  result 
in  much  more  significant  differences  in  degree  of  cure,  and,  thus,  in  instantaneous  modulus  aixl  percent 
shrinkage,  than  in  the  graphite/epoxy  system.  Consequently,  we  anticipate  a  greater  potential  for  stress 
developmem  in  the  thick  glass^lyester  lamiruttes  due  to  itt  sensitivity  to  thennally  induced  gradients  in 
degree  of  cure. 

3.3  Volumetric  Resin  Shrinkage  Effects  on  Stress  DeveloumenL 

3.3.1  Graphite/Etwxv.  Figure  3  shows  the  predicted  residual  transverse  stress  profiles  through  the 
thickness  of  a  17.78-cm  unidirectional  graphite/epoxy  laminate.  Profiles  are  shown  as  a  function  of  the 
magnitude  of  resin  chemical  shrinkage  during  cure.  Note  that  the  profile  for  the  0.0%  cure  shrinkage  case 


8 


Table  3.  Fiber  and  Matrix  Constituent  Properties 


Property 

Glass 

Polyester 

Graphite 

Epoxy 

E,(MPa] 

7.308  X  10^ 

eqn.  1 

eqa  1 

EifMPa) 

eqa  ! 

eqa  1 

Vu 

0.22 

0.40 

0.20 

0.33 

VjS 

0.22 

0.40 

0.20 

0.35 

0.22 

0.40 

0.30 

0.35 

0,2[MPa] 

2.992  X  10* 

eqn.  2 

2.738  X  10* 

eqa  2 

G.jCMPa) 

2.992  X  10* 

eqn.  2 

2.738  r.  10* 

eqa  2 

GafMPa] 

2.992  X  10* 

eqn.  2 

2.738  X  10* 

eqa  2 

Or 

3.04  X  KT* 

7.20  X  la* 

-9.00  X  KT^ 

5.76  X  10'^ 

02 

3.04  X  10^ 

7.20  X  la* 

7.20  X  10^ 

5.76  X  10-’ 

Table  4.  Resin  Characteristics  During  Cure 


1  Property 

Polyester 

Epoxy 

2.757 

3.447 

E„," 

2.757  X  10* 

3.447  X  10* 

v’’ 

''  «b 

4.0  -  6.0 

1.0  -  3.0 

does  not  represent  the  rninimum  residual  transverse  stress,  as  there  is  a  definite  contribution  from  thermal 
expansion  of  the  fiber  and  resin.  The  authors  propose  that  the  residual  stress  profile  induced  by  thermal 
expansion  effects  alone  is  of  the  same  magnitude  but  in  direct  competition  with  the  stresses  induced  by 
chemical  shrinkage  during  cure  for  the  graphiie/epoxy  system  being  examined. 

Figure  4  compares  the  residual  transverse  stress  profiles  induced  by  the  isolated  effects  of  thermal 
expansion  and  chemical  shrinkage  in  a  17.78-cm  laminate.  Thus,  profiles  represent  either  the  effects  of 
thermal  expansion  only  (CUE)  or  of  volumetric  shrinkage  only  (V ■  1.0,  3.0,  and  6.0%).  We  see  that 
the  magnitude  of  the  stresses  developed  by  the  thermal  expansion  strains  is  very  close  to  that  developed 
by  the  1.3%  volumetric  shrinkage  case.  As  the  chemical  cure  shrinkage  of  the  3501-6  epoxy  is  thought 
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Figure  2.  Normalized  Resin  Volumetfic  Shrinkage  During  Isothermal  Cure 
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3.3.2  Glass/Tolvester.  A  typical  value  for  the  cure  shrinkage  of  the  4102  polyester  resin  is 
approximately  6.0%.  as  opposed  to  l.S%  for  the  3S01'6  epoxy.  The  increased  contribution  of  the 
volumetric  shrinkage,  as  well  as  the  increased  severity  of  the  cure  kinetics,  leads  us  to  expect  significantly 
higher  residual  stress  development  in  the  glass/polyester  system  as  compared  to  the  graphite/epoxy 
system.  Figure  5  illustrates  process-induced  transverse  stress  profiles  for  a  254-cm  unidirectional  glass/ 
polyester  laminate  as  a  ftmetion  of  volumetric  shrinkage.  The  reversal  of  the  parabolic  stress  profiles  is 
a  result  of  the  interactions  of  thickness  and  cure  kinetics  inducing  a  predomirumtly  inside-to-outside  cure 
rather  than  the  outside-to-inside  cure  exhibited  by  the  tiiicker  graphite/epoxy  sections.  It  is  again 
demonstrated  that  the  case  with  0.0%  resin  shrinkage  does  not  minimize  process-induced  stress 
development.  Here  it  is  demonstrated  that  a  finite  resin  shrinkage,  on  the  order  of  1.0%,  produces  the 
nearly  stress-free  laminate.  Note  that  the  peak  predicted  residual  stresses  in  this  2.S4-cm  glass^Iyester 
laminate  *  6%)  are  significantly  larger  than  those  predicted  in  the  17.78-cm  graphite/epoxy  laminate 
=  1.5%).  This  demonstrates  the  importance  of  volumetric  resin  shrinkage  and  cure  kinetics  on  the 
process-induced  residual  stress  development  in  thick-section  laminates. 

Figure  6  shows  the  predicted  longitudirutl  residual  stress  profiles  as  a  function  of  volumetric  shrinkage 
for  the  same  laminate.  The  magnitude  of  the  residual  stress  in  the  fiber-dominated  longitudinal  direction 
is  significantly  lower  than  in  the  transverse  direction,  where  processing  strains  are  high.  Also,  it  is  worth 
noting  that  the  volumetric  shrinkage  here — which  most  nearly  produces  a  stress-free  laminate  in  the 
longitudinal  direction  (approximately  2.0%)— is  significantly  larger  than  that  for  the  transverse  direction 
(approximately  1.0%).  This  indicates  that  it  is  inherently  impossible  to  tailor  the  volumetric  shrinkage 
of  a  resin  to  produce  stress-free  profiles  in  both  the  trartsverse  and  longitudinal  directions.  As  such,  it  is 
more  desirable  to  tailor  the  resin  volumetric  shrinkage  to  minimize  the  residual  stress  development  in  the 
weak  transverse  direction,  where  matrix  cracking  could  more  easily  develop. 

Figure  7  compares  the  transverse  residual  stress  profiles  induced  by  the  isolated  effects  of  thermal 
expansion  strains  (CTE)  to  the  profiles  resulting  from  volumetric  shrinkage  only 
(V'^t  «  I.O,  3.0,  and  6.0%)  for  the  same  2.S4-cm  glass/polyester  laminate.  The  thermal  expansion  effect 
is  again  seen  to  be  opposite  in  sign  to  the  shrinkage  effects  and  of  similar  magnitude  to  the  1.0% 
shrinkage  case.  However,  as  the  polyester  resin  actually  undeigoes  approximately  6.0%  shrinkage  during 
cure,  the  effects  of  thermal  expansion  are  nearly  insignificant  with  respect  to  the  total  residual  stress 
development  in  this  example. 
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Figure  7.  Thermal  Expansion  and  Resin  Shrinlcagc  (Plotted  Scparaidv)  Coniribuiions  on  ihe  Predicted  Transverse  Residual  Stress 
Profiles  in  a  2 J4-an  GI«ss/PDlvesier  Lamiime. 


Predicted  residual  transverse  stress  profiles  for  a  S.08<m  glass^lyester  lamiiiaie  are  illustrated  as 
a  function  of  the  magnitude  of  volumetric  shrinkage  in  Figure  8.  Note  the  reversal  of  the  parabolic  stress 
profiles  from  the  2^-cm  laminate.  This  results  from  the  faa  that,  in  this  case,  the  cure  froiu  propagates 
from  the  outside  to  the  inside  of  the  laminate  due  to  the  interactions  between  the  tMdoiess  and  the  thermal 
history  of  the  laminate.  It  is  interesting  to  note  that  the  magnitude  of  the  residual  transverse  stresses  in 
the  5.08-cm  laminate  is  less  than  that  of  the  laminate.  This  is  probatdy  attributatde  to  the 

competing  effects  of  simultarvous  cure  from  propagations  (outside  to  inside  and  inside  to  outside),  as  is 
evidenced  by  the  comf^ex  stress  fields  through  the  laminate  thickness. 

It  is  interesting  to  note  that  the  magnitude  of  the  diemical  shrinkage  required  to  counteihalance  the 
residual  transverse  stresses  induced  by  thermal  expansion  effects  for  a  given  system  is  not  constant  and 
is  dependent  upon  the  comfdex  interactions  between  cure  kinetics,  laminate  thickrtess,  and  processing 
history.  From  the  inspection  of  Figure  S  it  is  appatent  that  for  a  2  J4<cm  laminate  (glass^nlyester).  a 
volumetric  shrinkage  of  slightly  over  1 .0%  will  minimize  transverse  residual  stress  distributions.  However, 
cbse  examination  of  Figure  7  reveals  that' for  a  S.08*cm  laminate  of  the  same  system,  a  shrinkage  of 
slightly  less  than  1.0%  is  necessary  to  minimize  the  transverse  residual  stresses. 

3.4  Experimental  Correlation.  Hahn  and  Kim  (1990)  subjected  several  grafrtiite/epoxy  (T30(V3S01-b) 
unidirectional  and  unsymmetrical  (Q^]  laminates  to  iruemipted  cure  cycles  to  monitor  the  development 
of  transverse  modulus  (in  the  unidirectional  laminates)  and  curvature  irxluoed  by  cure  shrinkage  and 
thermal  strains  (in  the  unsymmetrical  laminates).  The  laminates  were  approximately  1  mm  in  thickness 
and  the  cure  cycles  followed  the  manufacturers  recommended  cure  cycle  (MRC).  At  each  of  a  series  of 
points  during  the  MRC.  a  unidirectional  and  a  cross-ply  laminate  were  cooled  to  room  temperature  at 
3°  CVmin.  The  unidirectional  laminates  were  tensile  tested  to  determine  the  transverse  modulus,  artd  the 
deflections  of  the  unsymmetrical  laminates  were  measured  to  determine  process-induced  curvatures. 
Figure  9  shows  excellent  agreement  between  the  measured  and  predicted  transverse  modulus  developmeru 
based  on  the  material  model  presented  herein.  Note  that  the  horizontal  axis  refers  to  the  time  of  cute  stop. 
This  indicates  the  time  at  which  the  MRC  was  interrupted  and  cooling  to  room  temperature  initiated  and 
does  not  indicate  the  time  of  the  full  cycle  to  which  the  laminate  was  subjected,  as  not  all  laminates  were 
cooled  from  the  same  temperature. 

Figure  10  shows  the  development  of  nondimensionalized  curvature  (with  reqrect  to  laminate  thickness) 
as  a  function  of  the  time  of  cure  stop.  Again,  we  see  good  agreement  between  measured  and  predicted 
values.  Significant  over-prediedon  of  the  curvature  earty  in  the  cycle  is  thought  to  arise  from  the  highly 
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Hgluc  8.  Predicted  Transverse  Residual  Stress  ProrUes  «s  i  Function  of  Magnitude  of  Resin  Shrinkage  in  a 
fThemial  Expansion  and  Resin  Shrinkage  Effects  Superimposed). 


igure  9.  Experimemal  vs.  Prediaed  Composite  Transvcfse  Modulus  Devciopmcm  in  an  8-p1v  UnidirecUonal  Granhite/FjMKY  i  -iiininai<> 


viscoelastic  behavior  of  the  matrix  which  has  not  yet  reached  gelation.  Also,  the  fracture  surfaces  of  the 
tensile  tested  specimens  in  the  region  of  the  curvature  over  prediction  exhibited  little  or  no  interfacial 
bonding  between  the  fibers  and  the  matrix,  leading  to  a  relatively  unconstrained  matrix  and  reduced  stress 
development  By  defining  the  gel  point  as  the  point  at  which  the  viscosity  of  the  resin  reaches  100  Pa>s, 
Hahn  and  Kim  (1990)  note  that  the  curvature  and  the  transverse  modulus  both  increase  sharply 
immediately  after  the  gel  point  possibly  indicating  that  the  developmem  of  residual  stresses  can  be 
significanfiy  affected  by  the  mechanical  properties  of  the  resin  before  fiill  cure  has  been  reached. 

4.  CONCXUSIONS 

A  process-induced  stress  model,  employing  a  cure  simulation  analysis  and  laminated  plate  theory,  was 
used  to  study  the  evolution  of  macroscopic  in-plane  residual  stresses  during  the  cure  of  thick-section 
glass^lyester  and  graphite/epoxy  laminates.  Due  to  the  neglect  of  thermal  relaxation  effects,  model 
predictions  are  expected  to  provide  upper  bound  estimates  on  residual  stress  developmeiu.  Material 
modeb  were  proposed  to  describe  chemical  hardening  and  shrinkage  of  the  resin  as  a  function  of  degree 
of  cure.  The  effects  of  cure  kinetics  and  changes  in  specific  volume  of  the  resin  during  cure  and  cool 
down  were  found  to  be  significant  driving  forces  for  the  development  of  residual  stresses.  The  model  was 
employed  to  investigate  the  competing  mechanisms  of  cure  shrinkage  and  thermal  expansion  in  residual 
stress  development  in  thick-section  laminates.  It  was  demonstrated  that  the  simultarteous  minimisation 
of  residual  stress  in  both  the  longitudinal  and  transverse  directions  through  tailoring  of  the  resin  shrinkage 
is  not  possible.  Furthermore,  the  optimum  resin  shrinkage  to  minimize  residual  stresses  in  the  weak 
transverse  direction  is  strongly  dependent  on  reaction  kinetics,  specimen  thickness,  and  the  cure  cycle 
employed  during  processing.  Correlation  of  model  predictions  with  previously  published  experimental 
data  is  in  good  agreement  This  correlation  provides  confideiKe  in  the  material  models  proposed  to 
describe  the  resin  behavior  during  cure  and  adds  credibility  to  the  methodology  developed  for  predicting 
the  evolution  of  residual  stress  and  deformation  in  thick-section  therroosei  laminates  during  cure. 
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